Graphene: A Pseudochiral Fermi Liquid 
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Doped graphene sheets are pseudochiral two-dimensional Fermi liquids with 
abnormal electron-electron interaction physics. We address graphene's Fermi 
liquid properties quantitatively using a microscopic random-phase-approxima- 
tion theory and comment on the importance of using exchange-correlation po- 
tentials based on the properties of a chiral two-dimensional electron gas in 
density-functional-theory applications to graphene nanostructures. 
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The seminal work of Novoselov et al. (1, 2) and Zhang et al. (3), which demonstrated electrical 
contact to an isolated single-layer graphite flake, has heightened interest in the graphene two- 
dimensional electron system (2DES). Graphene's honeycomb lattice (Fig. 1) has two-atoms per 
unit cell and its 7r-valence band and 7r*-conduction band touch at two inequivalent points in 
the honeycomb lattice Brillouin-zone. The energy bands near K and K' are described at low 
energies by a spin-independent massless Dirac equation (4): 

H=VT (Oipi + (T2P2) , (1) 

where r = ±1 for K and K' valleys, the Pi are envelope function momentum operators, and the 
<Ji are Pauli matrices which act on the sublattice pseudospin degree-of-freedom. Like the bands 
of semiconductor heteroj unction 2DESs, graphene's bands can be described by a continuum 
model with a single parameter, in this case the Dirac velocity v instead of the band mass m b . 
Just as a semiconductor provides a realization of non-relativistic quantum mechanics with a 
material-dependent mass m b , graphene provides a realization of relativistic quantum mechanics 
with a material specific velocity of light. When Coulombic electron-electron interactions are 
included, doped graphene represents a new type of many-electron problem, distinct from both 
an ordinary 2DES and from quantum electrodynamics. The Hamiltonian in Eq. (OQ) differs from 
a Schrodinger equation in two crucial respects: i) its spectrum is not bounded from below and 
ii) its eigenstates have definite projection of pseudospin along the direction of momentum, i.e. 
definite pseudochirality, rather than definite pseudospin. We refer to the graphene 2DES as the 
chiral 2DES (C2DES). In this article we explain why the C2DES and the ordinary 2DES have 
distinctly different Fermi liquid properties. 

2DESs have been a fertile source (5) of surprising new physics for more than four decades, 
most notably in the presence of a strong magnetic field when they exhibit the integer (<5) and 
fractional (7) quantum Hall effects. Although the exploration of graphene is still at an early 
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stage, it is already clear (1-3, 8-10) that the strong field properties of a C2DES are different 
from and as rich as those of a semiconductor heteroj unction 2DES. In the absence of a field, 
ordinary 2DESs are normal Fermi liquids (11, 12) in which interactions alter the Fermi veloc- 
ity v — > v*, introduce marginally irrelevant effective interactions between quasiparticles on the 
circular Fermi surface, and diminish the fraction Z of the spectral weight in the one-particle 
Green's function associated with its quasiparticle peak. The Fermi liquid phenomenologies of 
a C2DES and an ordinary 2DES have the same structure, since both systems are isotropic and 
have a single circular Fermi surface as illustrated in Fig. 2. The strength of interaction effects 
in an ordinary 2DES increases with decreasing carrier density. At low densities, the quasiparti- 
cle weight Z is small, the velocity is suppressed, the charge compressibility changes sign from 
positive to negative, and the spin- susceptibility is strongly enhanced. These effects, described 
with reasonable consistency (13-19) by theory and experiment, emerge from an interplay be- 
tween exchange interactions and quantum fluctuations of charge and spin in the 2DES. In the 
C2DES we find that interaction effects also strengthen with decreasing density, although more 
slowly, that the quasiparticle weight Z tends to larger values, that the velocity is enhanced 
rather than suppressed, and that the influence of interactions on the compressibility and the 
spin-susceptibility changes sign. These qualitative differences are due to exchange interactions 
between electrons near the Fermi surface and electrons in the negative energy sea, to quasi- 
particle chirality, and to interband contributions to C2DES charge and spin fluctuations. The 
interband excitations are closely analogous to virtual particle-antiparticle excitations of a truly 
relativistic electron gas. 

The technical calculation (20) on which our conclusions are based is an evaluation of the 
electron self-energy E of the C2DES near the quasiparticle-pole. E describes the interaction 
of a single electron near the 2DES Fermi surface with all states inside the Fermi sea, and with 
virtual particle-hole and collective excitations of the entire Fermi sea, as illustrated in Fig. 2. 
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As we discuss more explicitly below, a direct expansion of electron self-energy in powers of 
the Coulomb interaction is never possible in a 2DES because of the long-range of the Coulomb 
interaction. Our results for the C2DES are based on the random phase approximation (RPA) in 
which the self-energy is expanded to first order in the dynamically screened Coulomb interac- 
tion W (setting K — 1): 
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P^J (27T) 

where s = + for electron-doped systems and s = — for hole-doped systems, j3 = 1/(/cbT), 



(2) 



W{q, zfl) = v q + v 2 qXpP (q, ifi) , (3) 

is the RPA density-density response function, 1S i ts non-interacting limit (21-24), and 
e(q, iCl) is the RPA dielectric function. For definiteness, we limit our discussion to an electron- 
doped system with positive chemical potential fi: the Fermi liquid properties at negative doping 
are identical because of the C2DES model's particle-hole symmetry (see also Fig. 2). 

In Eq. © tu n = (2n + 1)tt/(3 is a fermionic Matsubara frequency, the sum runs over all 
the bosonic Matsubara frequencies Vt m = 2mir/ ' (3 while in Eqs. © and ©, v q is the bare 
unscreened Coulomb interaction in 2D, v q = 2ne 2 /(eq) where e is an effective dielectric con- 
stant. The first and second terms in Eq. © are responsible respectively for the exchange in- 
teraction with the occupied Fermi sea (including the negative energy component), and for the 
interaction with particle-hole and collective virtual fluctuations. The factor in square brackets 
in Eq. ©, which depends on the angle #k,k+ q between k and k + q, captures the dependence 
of Coulomb scattering on the relative chirality ss' of the interacting electrons. The Green's 
function G^(k, iuS) = l/[iu — £ s (k)] describes the free propogation of states with wavevec- 
tor k, Dirac energy £ s (k) = svk — /i (relative to the chemical potential) and chirality s = ±. 
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After continuation from imaginary to real frequencies, iuj — > uj + irj, the quasi-particle weight 
factor Z and the renormalized Fermi velocity can be expressed (20) in terms of the wavevector 
and frequency derivatives of the retarted self-energy E^ t (k, a;) evaluated at the Fermi surface 
(k = k F ) and at the quasiparticle pole uj = £+(k): 
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Following some standard manipulations (20) the self-energy can be expressed in a form conve- 
nient for numerical evaluation, as the sum of a contribution from the interaction of quasiparti- 
cles at the Fermi energy, the residue contribution S rcs , and a contribution from interactions with 
quasiparticles far from the Fermi energy and via both exchange and virtual fluctuations, the line 
contribution E lme . In the zero-temperature limit 
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E^(k,a;) = -5: 
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Note that at the Fermi energy <9 fc E^ s (k, uj) vanishes, and d w E+ s (k, uj) involves an integral over 
interactions on the Fermi surface that are statically screened. These expressions differ from the 
corresponding 2DES expressions because of the relative chirality dependence of the Coulomb 
matrix elements, because of the linear dispersion of the bare quasiparticle energies, and most 
importantly because of the fast short-wavelength density fluctuations produced by the interband 
contribution to x^HQi illustrated in Fig. 3. 
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Our results for Z and v*/v are summarized in Fig. 4 as a function of the C2DES dimension- 
less coupling constant (restoring h) 

2vre 2 e 2 
ek F evn 

The appropriate value of / for a particular graphene sheet is dependent on its dielectric en- 
vironment; for graphene on Si0 2 / ~ 2. As we discuss at greater length below, graphene's 
Fermi liquid properties depend only weakly on the carrier density which is expressed in these 
figures in terms of the cut-off parameter A. The trends exhibited in Fig. 4 can be understood by 
considering the limits of small / and the limit of large q at all values of /. In the former limit 
screening is weak except at extremely small q. In <9 w E^ s (k, uo), for example, the integral over 
q diverges logarithmically at small q when e(q, uo = 0) is set equal to one, i.e. when screening 
is neglected. Screening cuts off this logarithmic divergence at a wavevector proportional to / 
so that c^E+^k, u) has a contribution proportional to / ln(/) at small /. Because e(q, uo = 0) 
happens to be independent of q for transitions between Fermi surface points, it is possible to 
evaluate <9 w E^ s (k, u) analytically. We find that 
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Similar small q f ln(/) contributions appear in the other elements which contribute to Z and 
v* . All this behavior is very familiar from the case of the normal 2DES; the new differences 
present in the chiral C2DEG are ones of detail. At large q, on the other hand, interband charge 
fluctuations dominate e(q,u;) — 1, which approaches its simple undoped system form. It be- 
comes especially clear when uo is expressed in units of vq that the typical value of e(q, uo) at 
large q is ~ 1 with a non-trivial dependence on /. The q integrals all vary as q' 1 , requiring 
that the C2DES model be accompanied by an ultraviolet cut-off which for the case of graphene 
should be (24) q c ~ l/a where a is the graphene lattice constant. Since the crossover between 
intraband and interband screening occurs for q ~ fc F , it follows that both <9fcE lme and ct,E lme 



have contributions that are analytic in / and vary as ln(A) where A = q c /k-p when A is large. 
To leading order in ln(A) we find that 
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Note that A(/) = //4-3vr 2 / 2 /256 + ... and $(/) = 1/3 -3vr 2 //256 + ... are analytic functions 
of / because interband polarization screening does not essentially alter the Coulomb interaction 
at large q. 

The asymptotic expressions (PTTj) and (fT2~l) approximately capture the contribution to the 
corresponding Fermi liquid parameters from interactions over the wavevector range from ~ 
kp to ~ q c . As the density decreases and kp — » this contribution dominates. The Fermi 
wavelength then acts like a cut-off on the renormalization group flows which appear in the 
theory (26) of interaction effects in undoped graphene. The fact that the velocity increases in 
this regime can be understood qualitatively using Hartree-Fock theory (24), which is accurate 
at small / when A is large. In Hartree-Fock theory the enhanced velocity is due to the reduced 
exchange energy in a right-handed band when the negative energy sea is left-handed. In Fig. 4 
we have also shown cut-off and coupling constant dependence of the antisymmetric Landau 
parameter F a , which is defined in terms for the antisymmetric Landau interaction function 
/ a (cos(p)) as (20) 

2 ^/ a (cos(^))cos(^), (15) 
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where is* = gk-p/{2nv*) is the density-of-states of the interacting system at the Fermi surface 
[/ a (cos(</?)) is obtained from the Fermi surface dependence of the self-energy (27)]. Physically, 
F a determines the spin susceptibility xs — (v/v*)/(l + F a ). From our results in panels (b) and 
(c) of Fig. 4 we predict a rather large suppression of the spin susceptibility which could be mea- 
sured in weak-field Shubnikov-de Haas magnetotransport experiments using a tilted magnetic 
field to distinguish spin and orbital response (14). 

Our findings have important implications for density-functional-theory (DFT) and tight- 
binding modeling of ribbons (28-30), quantum dots (31, 32), and other nanostructures made 
from graphene. Because of the pseudo-chiral properties of bulk quasiparticles, states tend to 
have a lower energy when they have the majority chirality. This interaction effect depends 
specifically on intersite coherence and is completely missing in the local-density-approximation 
and in other approximations for exchange and correlation potentials commonly used in DFT. 
The accuracy of graphene nanostructure electronic structure calculations would be improved if 
they used exchange-correlation potentials based on the properties of the C2DES rather than on 
the properties of the ordinary 2DES. 

This work has been supported by the Welch Foundation, by the Natural Sciences and En- 
gineering Research Council of Canada, by the Department of Energy under grant DE-FG03- 
02ER45958, and by the National Science Foundation under grant DMR-0606489. 
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Figure 1 : Honeycomb lattice of a single layer graphite flake with one sublattice in yellow and 
the other sublattice in blue. In the continuum limit the sublattice degree of freedom may be 
regarded as a pseudospin. When momentum k is measured away from the Dirac points at 
the K and K' Brillouin zone corners, band eigenstates have definite projection of pseudospin 
in the k direction, i.e. definite pseudochirality. The angle 0k above denotes the momentum- 
dependent phase difference between wavefunction amplitudes on the two sublattices. For spin- 
1/2 quantum particles this angle is the azimuthal orientation of a pseudospin coherent state in 
the equatorial plane. 
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Figure 2: In a weakly doped material, graphene's energy bands can be described by a massless 
Dirac equation in which the role of spin is played by pseudospin. Like an ordinary 2DES, 
doped graphene has a circular Fermi surface. The Fermi liquid properties of graphene are a 
consequence of both exchange interactions between quasiparticles near the Fermi surface and 
states in the positive and negative energy Fermi seas and of interactions with both intra-band 
(short red vertical arrow) and inter-band (long red vertical arrow) virtual fluctuations of the 
electronic system. The yellow arrows in this figure indicate the pseudospin chirality of band 
eigenstates. Because of the difference in chirality between positive and negative energy bands, 
the velocity of graphene quasiparticles is enhanced by inter-band exchange interactions, tending 
to protect the system from magnetic and other instabilities, and reducing both charge and spin 
response functions. 
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Figure 3: "Lindhard" function x^HQi ^) °f a C2DES, in units of the non-interacting density- 
of- states at the Fermi surface v = gk F /(2%v), as a function of q/k F and 0//ion the imaginary 
frequency axis. k F = (Ann/g) 1 ^ 2 is the Fermi wavenumber, [i = vk F the Fermi energy, n the 
electron density and the flavor multiplicity g = g s g v = 4 for graphene because of its two-fold 
valley degeneracy. Because of interband fluctuations x diverges linearly with q for q — > oo 
and decays only like for — > oo in the C2DES, in contrast to the q~ 2 and fi~ 2 behaviors 
of the ordinary 2DES. In the static = limit x (q> 0) = — v for all q < 2k F for both chiral 
and ordinary 2DESs. 
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Figure 4: Density and coupling constant / dependence of some C2DES Fermi-liquid param- 
eters. The density is specified by A = q c /k-p. The density range studied most extensively in 
experiment, n ~ 10 11 cm -2 to n ~ 10 13 cm" 2 , corresponds to A = 100 to A = 10. In all 
panels the black solid line corresponds to the highest value of the cut-off parameter we have 
considered, A = 2.7 x 10 5 . The red dashed line illustrates the RPA Fermi-liquid parameters of 
an ordinary non-chiral 2DES with parabolic bands. In this case the / = \f2 r s [see Eq. ©], 
where r s = (7rna|)~ 1//2 is the usual Wigner-Seitz density parameter and ae = eh 2 / (m h e 2 ) the 
effective Bohr radius. From the left the three panels show: (a) the quasiparticle renormalization 
factor Z evaluated from Eq. ©; (b) the velocity renormalization factor evaluated from Eq. ©; 
and (c) the £ = dimensionless Landau parameter Fq which characterizes spin-dependent 
quasiparticle interactions. The color coding for A is the same in all panels. 
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